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Separated flows past complex geometries are modeUed by discrete vortex techniques. The flows are 
assumed to be rotational and inviscid, and a new technique is described to determine the streamfunctions 
for linear shear proflles. The geometries considered are the snow cornice and the backward-facing step, 
' whose edges allow for the separation of the flow and reattachment downstream of the recirculation 
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, regions. A point vortex has been added to the flows in order to constrain the separation points to 



o 



be located at the edges, while the conformal mappings have been modifled in order to smooth the 
sharp edges and let the separation points be free to oscillate around the points of maximum curvature. 
Unsteadiness is imposed on the flow by perturbing the vortex location, either by displacing the vortex 
from equilibrium, or by imposing a random perturbation with zero mean on the vortex in equilibrium. 
The trajectories of passive scalars continuously released upwind of the separation point and trapped by 
the recirculating bubble are numerically integrated, and concentration time series are calculated at fixed 



o 

^ ' locations downwind of the reattachment points. This model proves to be capable of reproducing the 



trapping and intermittent release of scalars, in agreement with the simulation of the flow past a snow 
cornice performed by a discrete multi-vortex model, as well as with direct numerical simulations of the 
flow past a backward-facing step. The simulation results indicate that for flows undergoing separation 
and reattachment the unsteadiness of the recirculating bubble is the main mechanism responsible for 
the intense large-scale concentration fluctuations downstream. 

1 Introduction 

Transport and dispersion of passive scalars in turbulent boundary layers over complex geometries is an 
important problem of practical relevance, which has been widely investigated in recent years mainly 
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because of a growing concern about the environmental issue, but also because of its engineering ap- 
plications such as optimization of combustors or evaluation of mixing properties of flows in chemical 
reactors. 

However, the complexity of flows in which separation and reattachment occur make them less suitable 
to analytical studies than regular boundary layer flows. Relatively little is known about the behaviour 
of a phune near hills, or buildings, or, in general, over irregular terrain (Hosker 1984). For a classical 
problem like the flow over a backward-facing step, there is not an applicable mathematical technique 
for modeling the evolution of the scalar field. Laboratory experiments generally provide an Eulerian 
portrait of this flow (e.g. Bradshaw & Wong 1972), but with no insight into the dispersion process. More 
recently, Le, Moin & Kim (1997) performed a direct numerical simiilation of a backward-facing step for 
Reynolds numbers as high as 5100 revealing, among the others, some of the unsteady flow characteristics. 
They detected the oscillation of the (spanwise-averaged) reattachment location, in agreement with the 
results of several other works (e.g. Driver, Seegmiller & Marvin 1983, 1987). The shear layer composed 
of many small high-intensity vortices, which extends to the reattachment point, rolls up forming a 
large-scale structure in the recirculation region: the periodic detachment of this large-scale vortex from 
the step causes the oscillation of the reattachment point. The main aim of this paper is to model the 
unsteadiness of reattached flows, and to assess how it affects the dispersion of passive scalars released 
by a continuous point source. 

We will model separated flows by means of point vortex techniques, assuming the flow as two- 
dimensional and inviscid. We will consider two different geometries: a backward-facing step, and a 
"snow cornice" . Snow cornices are natural devices that control the flow separation on momitain crests 
by trapping vortices. The modeling technique is not new: Ringleb (1961) studied the steady separated 
flow past a snow cornice assuming a two dimensional potential flow, which allowed him to solve the 
problem analytically by means of the classical method of conformal mapping. To represent a snow 
cornice, Ringleb devised mappings to transform the real axis in the complex ^-plane onto a line forming 
a sharp edge in the complex z-plane. In this way the complex potential of a uniform flow over the real 
axis in the ^-plane was ascribed to the flow on the corresponding region in the z-plane. A point vortex 
(together with its reflected image) was added to the flow to model the recirculating region. The vortex 
was in equilibrium and satisfied a steady Kutta condition. In a more recent paper, Cortelezzi, Leonard 
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& Doyle (1994) apply point vortex techniques to model the unsteady separated flow over a semi-infinite 
plate, and to assess active control techniques. 

Because of the incompressibility condition, the Lagrangian motion of a passive scalar is governed by 
a Hamiltonian system whose conjugate variables are the Cartesian coordinates of the scalar, and whose 
Hamilton function is the streamfunction. The motion of a point vortex is governed by a Hamiltonian 
system as well, where the Hamilton function is related to the Green function for the Laplacian of the 
flow region (Masotti 1931; Lin 1941). 

In this respect, the motion of the vortex in the Ringleb model is integrable since it is governed by a 
one-degree of freedom autonomous Hamiltonian system. Incidentally, although we found that Ringleb 
made a formal error in determining the law of the motion of a vortex in the presence of a wall, and 
therefore failed to derive the correct equilibrium conditions, his results are qualitatively correct (see 
Appendix B). 

The streamfunction, namely the Hamiltonian of the fluid particles, depends on the vortex position 
and therefore is in general time-dependent (non-autonomous). For such a system we expect a non 
integrable, chaotic, particle motion (Novikov & Sedov 1979; Aref 1983). However, if the vortex is 
located at its equilibrium point, the Hamiltonian of the particles is autonomous and the motion is 
integrable and regular. For the integrable case, the flow field is virtually divided in two main bodies: 
the fluid entrained by the vortex (closed streamlines encircling the vortex), and the free-stream flow 
(open streamlines). The two regions do not exchange fluid with each other and are separated by a 
streamline which leaves the sharp edge of the wall, and reattaches downstream (i.e. a heteroclinic 
orbit). For instance, in the transformed ^-plane this flow is a vortex pair, where the real axis is the line 
of symmetry. 

The present model is different from Ringleb's in three main aspects: 

(i) We assume a shear flow as the asymptotic upstream condition in order to have a more meaningful 
flow than the potential one, but still simple enough to allow for an analytical study. 

(ii) The system has been made unsteady by perturbing the equilibrium of the vortex. 

(iii) Finally, we consider geometries where the sharp corners are smoothed to blunted edges with 
large but flnite curvature, to satisfy an unsteady Kutta condition for a vortex of constant circulation 
not in equilibrium. 
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The large curvature edge still causes separation, in that according to the Kutta condition it is 
possible to define a vortex whose strength produces separation. A vortex slightly displaced from its 
stable equilibrium moves along a small periodic orbit, and caiises the oscillation of the detachment 
point around the maximum curvature point. Since the edge is smooth, the flow around the edge does 
not have any singularities and the Kutta condition can be considered fulfilled. 

A typical representation of the mean flow vertical profile by power laws is used in the atmospheric 
boundary layer in stable conditions. For instance, a linear shear profile has been considered by Hunt, 
Leibovich & Richards (1988) to study how the shear affects the flow past hills with low slopes. 

On the transformed (^-plane, the flow reduces to the oscillating vortex pair (OVP). Rom-Kedar, 
Leonard & Wiggins (1990) studied the OVP flow extensively, and elucidated the mechanisms of the 
phase-space Lagrangian chaos. The stable and unstable manifolds of the heteroclinic orbit intersect 
each other, and drive the fluid to be entrained and detrained by the vortex. The Eulerian counterpart 
of this phenomenon is the chaotic mixing of the fluid in the region surrounding the separatrix streamline. 
We are interested in the effects of these phenomena on the dispersion of scalars. 

We investigated how the separation and reattachment of an unsteady flow influence the behaviour 
of effluent plumes released near or inside the recirculating region. We found that the presence of the 
recirculation bubble dramatically affects the large-scale dynamics of passive scalars, causing important 
fluctuations of the concentration field downstream. 

In the following sections the Ringleb model for steady potential flow is recalled, then a method is 
proposed to study analytically unsteady rotational flows with constant vorticity. Dynamical systems 
formalism is used to flnd equilibrium configurations, to discuss their stability and to recognize and 
quantify chaotic mixing. Numerical simulations for two different geometries are presented, along with 
the concentration time series at several sampling points. 

2 Mathematical formulation of the steady potential flow 

In this section we describe a mathematical model of recirculating flow over complex geometries such as 
a snow cornice and a backward-facing step. This model, which was proposed by Ringleb (1961), is able 
to represent a two-dimensional steady separation from the edge. It constitutes the basis for the model 



4 



of unsteady, recirculating shear flow, which we present in §3. 

We consider a region bounded by a piecewise analytical curve in the z-plane, where z — x + iy. 
The solid boundary extends to infinity along the x-axis, and is characterized by a sharp corner. As 
a consequence of the Riemann mapping theorem (see, e.g., Nehari 1975; Henrici 1974), such a region 
can be obtained as the conformal image of the half-plane > in the ^-plane, where C = ^ + irj. The 
infinity in the z-plane corresponds to the infinity in the C-plane. For instance, the function suggested 
by Ringleb (1961), 

= = C + ^ (1) 

where Ci = + ir]i is a complex constant with 771 < 0, yields to the geometry plotted in figure la. The 
backward-facing step shown in figure lb can be obtained by the Schwartz-Christoffel mapping: 



^Vc^+iog c + Vc^ • (2) 



The model consists of a vortex of constant strength which satisfies a Kutta condition, in equilibrium 
with a steady free-stream velocity. The complex potential of the fiow is built by superimposition of 
basic flows: in the C-plane it is expressed by the function: 

w = q^C+^log(^) (3) 



where q^o is the free-stream velocity, and the second term on the right hand side is the complex potential 
of a point vortex of strength 7 located at Co and of its image at the complex conjugate position Q. 
Equation (^) represents the complex potential of a vortex pair of opposite signs in the (^-plane (see figure 
7). 

Since the complex potential is invariant under a conformal mapping z = z{(), the constant qoo is 
equal to the free-stream velocity Qoo in the z-plane: 

dw dw ( dz\^^ 

Qoo = 1™ -1- = 1™ -77 ^77 = 9oo ■ (4) 

z^oo dz (^00 dQ \dQ J 

The trajectories of a vortex in the physical z-plane are the solutions of the Hamiltonian system: 

dH . dH 

= yo = (5) 

oyo 0x0 
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where the Hamihonian H can be derived from the Hamiltonian H' of a vortex in the ^-plane, according 
to the Routh rule (e.g. Clements 1973): 



H = H' + — log 
An 



(6) 



where: 



-ff' = goo?7o + :i-log??o ■ (7) 

47r 



Derivations of the Hamiltonian H' and the Routh transformation rule are given in Appendix A. 
2.1 Steady Kutta condition 

The separation of a flow around a sharp or a large curvature corner is captured by the model if it 
satisfies the Kutta condition, that is, a relationship between the strength of the vortex and its position 
in order to eliminate the singularity of the flow. Therefore, the condition of zero velocity at the edge is 
imposed: 

where Cc indicates the location of the corner in the C-plane. 

For the complex potential w{C) @), the Kutta condition is satisfied when 

goo '70 

where we notice that there are infinite possible values for the ratio 7/(700- For the mapping = 0, 

whereas = — 1 for (||). 

On the physical plane, the complex velocity at the corner Zc is 

lim — — = lim —r--^ (10) 

z^zc dz C-»Cc aC dz 

whose value is finite when the corner is a cusp, as in the case of mapping (|l|), while it is null in the other 
case. 

In order to complete the modeling of the recirculating flow, we are now concerned with finding 
the equilibrium conditions for the vortex. We recall that the phase space of the dynamical systems 
(^ coincides with the flow region, in the sense that the Cartesian coordinates of the vortex (xq, yo) 
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represent the conjugate variables, and the level lines of the Hamiltonian H form the pattern of the 
possible trajectories of a vortex of assigned strength. Therefore the equilibrium locations (io = 0, 
yo = 0) are the fixed points of this map: the elliptic points correspond to stable equilibrium, the 
hyperbolic ones to unstable equilibrium. The model of a steady recirculating flow behind a corner is 
represented by a vortex in stable equilibrium (i.e. a trapped vortex) ^ and satisfying the Kutta condition. 
The solutions to the equilibrium equations, if any, depend on the geometry of the solid boundary and 
on the ratio of the vortex intensity to the asymptotic velocity, i.e. 7/^00 ■ Since the geometries we 
considered were found to allow vortex capturing, we expect a set of stable solutions. 

Therefore, the vortex coordinates (i^o, ??o) on the C-plane have to satisfy simultaneously the equilib- 
rium condition and the Kutta condition (^). They are the solutions of the equation: 



2.2 Comparison with a multiple vortex method 

The accuracy of the above vortex model has been checked by comparing some of the results with 
those generated by a more sophisticated vortex method. We used a discrete vortex method (Ceschini 
1993; Ferlauto 1996) to describe the impulsive start of the flow characterized by vortex shedding. We 
performed the simulation for the snow-cornice geometry. The method consists of assuming the initial 
flow around the wall as potential, and satisfying an unsteady Kutta condition by adding a free point 
vortex of suitable strength in a prescribed location close to the corner. The vortex moves because of, 
and according to, the wall and the background flow. The motion is then numerically integrated. At 
given time intervals a new vortex is introduced into the flow to restore the regularity. Figure 2 shows 
the pattern of the vortices shed throughout the transient, where the vortices strengths are proportional 
to the symbol sizes. The transient is characterized by a vortex sheet rolling up around the point of 
highest absolute circulation. This result is consistent with the description of the transient as given by 
Cortelezzi, Leonard & Doyle (1994) and Le, Moin & Kim (1997) for other vortex trapping flows. The 
vortex sheet increases in size and circulation until eventually reaching a more steady configuration, with 
small oscillation of its global circulation and centroid. 




(11) 



The vortex strength is then obtained by the Kutta condition (H). 
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We compared the asymptotic circulation and the centroid location with the strength and location 
of the single point vortex used in our model (the centroid location is evaluated assuming the strengths 
as masses). We plot in figure 3a the total circulation as a function of time (solid line), as well as the 
circulation for the single vortex model (dotted line) . In figure 3b we plot the distance of the centroid for 
the multiple vortex method from the equilibrium point of a single vortex. The two models are in good 
agreement with each other as far as the intensity of the vorticity field and the location of the centre of 
vorticity are concerned, showing that the large-scale characteristics of the flow arc well captured by a 
single vortex model. Of course, the strength of the single vortex model is in its simplicity, which permits 
the analytical description of a large-scale mechanism of dispersion, as discussed in §4. 



3 Shear flow 

In this section we describe a method to model the separation and recirculation of a shear flow over 
complex geometries. 

We assume a rotational flow field, with vorticity V x q = w = const and with the following velocity 
profile in the far field: 

lim q ^ Uoa - Luy . (12) 
The motion of passive scalars is governed by the Hamiltonian dynamical system: 

. d^p . dip 

where the streamfunction ■0 represents the Hamilton function and the scalar coordinates (x, y) are 
canonical variables. The streamfunction ?/' is related to the vorticity lj by the linear Poisson equation: 

= . (14) 



To solve the Poisson equation (14) we developed an analytical procedure based on conformal mapping, 
which is inspired by the method described by Tsien (1943). 
Let us assume 

ip = i>Lj + ipp (15) 
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where the streamfunction is relevant to the shear flow u = —ujy, i.e. 

V'a;=-^W2/^ (16) 

and the function tpp is harmonic in the flow region. 
Since 

VVc. = -co (17) 

and 

VVp = (18) 

equation ( [T^ ) is satisfied. However, the problem has been reduced to the determination of the harmonic 
streamfunction ■0p which ensures the fulfillment of the boundary conditions, that is, the far field velocity 
equal to {uoc ~ coy) and the impermeability of the wall. 

Therefore, we seek a complex potential Wp whose imaginary part is ipp. We define the flow regions 
in the z-plane as the conformal representation of the half-plane 9(C) > according to the mappings: 

^ = C+ ,Jl , (19) 

and 

z = V(C + *^)' - 1 + log (c + ^'5 + V(C + - l) (20) 

which are the modified Ringleb mapping (^ and the Schwartz-Cristoffel mapping (|^), respectively, 
where we introduced the real positive parameter 5 in order to smooth the corners (the smaller 5, the 
sharper the corner). 

The further mapping: 

transforms the half- plane 3(C) > onto the interior of the unit circle |A| = 1 in the A-plane. Summing 
up, the chain mapping z Q ^ \ maps the physical boundary of the flow region onto the unit circle in 
the A-plane, where A = — 1 corresponds to infinity in the z-plane. 



On the physical plane, the shear fiow (16) induces the following velocity component v normal to the 
boundary: 

V — uj yb sin/3 (22) 
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where the subscript b denotes boundary, and f3 — arg {dz / dC) is the angle between the boundary and 
the X-axis. 

According to equation (p2[), the condition of inipermeabihty for the boundary in the z-plane is 
satisfied when the velocity component normal to the boundary, due to the complex potential Wp, is 
equal to —v, that is: 



dw 



dz 



which translates on the A-plane as 







dz 








d\ 


|A| = 1 



(23) 



(24) 



As there are no point vortices (they can be added to the flow later), the complex velocity dwp/dz 
cannot be singular inside the flow fleld. Since the mappings z — > A are regular inside the unit circle 
l-^l 5: 1 (-^ 7^ ^1) in the A-plane, then the complex function dwp/dX must not be singular inside the unit 
circle as well, except for the point A = — 1 that corresponds to z = cxd. Therefore we can assume: 



dw. 
'dX 



^ =F(A) + ^(a„-z6„)A"-i 



(25) 



n=l 



where F{X) contains the possible singularities located in A = —1, that is, in general, 

oo 

F(A) = ^c,(A + l)-^". 



(26) 



To determine the Cj coefflcients, we recall that the velocity at infinity (z = oo), due to the complex 
potential Wp, is Uoo- Since for the mappings (|l^) or (^0|) the derivative at infinity is 



lim — ^ = 1 

z^oo dz 



we have 



lim — - = lim 

C^oo d(, 



dwp dX 
~dXdC 



lim 



dwp (A + 1)^ 
"dA 



-2i 



(27) 



(28) 



which implies that all the coefficients c,- for j > 2 are null, and that 



C2 



-2i Uc 



(29) 
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We determine the coefficient ci by considering that the mass flow Q across the boundary, due to 
the shear flow, has to be balanced by the potential flow. We have 

^+00 



r^ f '^y^ ^ ^ r 21 



21 +°° 



and also, according to the invariance property of the complex potential. 



The function -F'(A) is then fully determined: 

F{\) = -2z 



-.1 dX 



1 



TT Cl . 



Q 1 



(30) 



(31) 



(32) 



(A + 1)2 TT A+1 

According to equation (^0|), the mass flow Q relevant to the mapping (|9|) is null, while for the mapping 
(1^) it is 

Q--|(7r2 + 27r S) . 



equations (24) and (pa) evaluated at |A| — 1, that is A = exp(ii9): 



To calculate the second term on the right hand side of equation (^5|) , we combine equation ( psj ) with 

Q A 



(a„ cos ni5 + 5„ sin ni?) = —v 



n=l 



dz 
dX 

dz 
dX 



2tt 



X 



TT A+1 



(33) 



Since the right hand side of ( |33| ) is not singular, the coefficients a„ and 6„ on the left hand side tend to 
zero for large n. The Fourier series on the left hand side can be suitably truncated and the coefficients 
a„ and 6„ evaluated numerically. The FFT algorithms are efficient tools for this purpose. 
The complex potential Wp is finally obtained by integrating equation (p5[): 



2i u, 



A + 



^-^log(A + l) + 5: 

1 TT 



(34) 



In general, for geometries defined by the mappings (19) and (pO[), the fiow does not satisfy the Kutta 
condition, i.e. it does not separate at the corner. The rotational flow field does exhibit a recirculating 
bubble, but the separation point is unphysically located downstream from the high curvature corner. 
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To enforce the Kutta condition, a vortex can be added to the flow, as shown in §2 for a potential flow: 
the same considerations regarding its strength and equilibrium hold for a shear flow. 
By adding a point vortex to the flow, the complex potential becomes 



w„ 



2z W-oc 

A + f 



log(A + l) + ^ 



ib„ 



n=l 



7 

2i TT 



log 



C-Co 



(35) 



The complete streamfunction, relevant to the sum of the shear and the potential component of the 
flow, is 



which provides the following flow velocity: 

oo 

F(A)+^(a„-i6„)A"-i 



z = -ojy- 



The vortex velocity can be obtained as 



dX 
dz 



1 



dz 



io = -cj yo + hm 

Z~*ZQ 



dWr, 



7 



1 



dz 2TTi z — zo ^ 

and is governed by a Hamiltonian system as in the case of the potential flow discussed in §2: 



Xo 



dH 
dyo 



2/0 



where the Hamilton function H is: 



2 ^ 



2i Q 

— — log A + 1 

A+f TT 



E 



dH 

dxo 



A" 
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log ?7o - log 



(36) 



(37) 



(38) 



(39) 



(40) 



The strength and location of the vortex is then determined according to equations ( |38D and (|3£ 
imposing the flow at rest on the corner Zc and the vortex at its equilibrium location (io — yo = 0). 
Figure 4 shows streamlines and velocity vectors for a separated shear flow for a geometry given by the 
map ([fOl); the flow field is characterized by uj/uoo ~ —10 and j/uoo — —2.243. 

The stability of the equilibrium can be inferred from figure 5, which displays the contour levels of 
the Hamiltonian of the vortex, namely the possible trajectories of the vortex. The equilibrium point, 
marked with a square, coincide with an elliptic, i.e. stable, fixed point in the phase space. 
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3.1 The unsteady flow 

We examine the Lagrangian transport of passive scalars in the case of unsteady flow. Recalhng equations 
( p^ and (|3^), the motion is governed by the system of equations: 

d d 
X = ^'4!{x,y;xo,yo); y = -—ij{x,y;xo,yo) , (41) 

where the streamfunction ip depends on both the particle coordinates {x,y) and the vortex location 
{xq, yo). For a vortex located at its equilibrium position {xq , y(f ), the one-degree of freedom Hamiltonian 
system (^ij) is autonomous, and hence integrable. In this case, the phase portrait coincides with the 
streamline pattern for the steady flow as shown in figure 4. The phase space is characterized by a 
streamline connecting two stagnation points, or according to dynamical system jargon, a heteroclinic 
orbit connecting two hyperbolic fixed points. Figure 6a shows this streamline in the physical z-plane. 
Such a feature has great relevance because it is capable of triggering chaotic mixing if the system 
undergoes a time-dependent perturbation, as shown below. 

If the vortex is displaced from its equilibrium location {xq , y^) by a small quantity e, its trajectory 
follows a closed orbit, as shown in figure 5. Its periodic motion obeys equations of the form: 

xo = x^ + f{xl^,y^,s,t), yo^y^ +g{x^,y^,e,t) . (42) 

In this case, system ( p| ) is no longer autonomous, because now the streamfunction ( |36| ) depends on 
time through the time-dependent vortex coordinates (a;o(t), 2/o(0)- Formally, we can write 

■0 = ■4j{x,y;xQ,y^,e) 

- *^<---» SI) +°<^^)' («' 

and the governing equations (^) reduce to the time-dependent perturbation of an integrable Hamilto- 
nian system: 



(44) 



dy dy \dxo de dyo de 
d'tjj^ d / dip df dip dg 
^ dx dx \dxo de dyo de 

As shown in figure 6a, for e = the unstable manifold leaving the hyperbolic fixed point P~ 
joins smoothly the stable manifold going into the other fixed point , and the heteroclinic orbit of 
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the unperturbed system is a single line connecting and P+. However, it is well known that e ^ 
dramatically affects the phase portrait in that the time dependence of the system, if not associated to any 
kind of symmetry, causes non-integrability (Novikov & Sedov 1979; Aref 1983). The stable and unstable 
manifolds do not join smoothly but intersect each other transversally forming a paradigmatically chaotic 
tangle (e.g. Guckcnheimcr & Holmes 1983, Tabor 1989, Wiggins 1992). Since any intersection point 
maps into another intersection point, the number of transverse intersections is infinite, whereas the 
distance between intersection points tends to zero as they approach a hyperbolic fixed point. At the 
same time, the system preserves the area of the lobes embraced by the arcs of the manifolds clipped 
between two intersections. Figure 6b shows a Poincare section of the flow over the snow cornice, 
enlightening the scenario of the stable and unstable manifold intersecting transversally. 

The chaotic character of a system can be deduced whenever evidence of such transverse crossings is 
provided. For instance, several authors used the Melnikov technique (Melnikov 1963) to detect analyt- 
ically the transverse crossings between stable and unstable manifolds for (time-periodically) perturbed 
integrable systems (e.g. Rom-Kedar, Leonard & Wiggins 1990, Zannetti & Franzese 1994, Del Castillo- 
Negrete 1998 and references therein). 

4 Results 

The representation on the transformed ^-plane of a steady shearless flow over a step or a snow cornice is 
similar in many respects to a vortex pair flow, except for the differences due to the Routh correction (^), 
which affects the vortex velocities. The streamline pattern of a vortex pair flow in the frame moving 
with the vortices is shown in figure 7, where we can see the two hyperbolic fixed points (P^,P+) 
connected by heteroclinic orbits. Rom-Kedar, Leonard & Wiggins (1990) perturbed the vortex pair 
by an external periodic strain-rate field, and studied extensively the chaotic dynamics of the resulting 
Oscillating Vortex Pair flow (OVP). They focused on the dynamics of the lobes enclosed by the perturbed 
heteroclinic orbits. Rom-Kedar, Leonard & Wiggins (1990), and Wiggins (1992), provided valuable tools 
for a quantitative evaluation of the mixing process in terms of the amount of fluid involved and the 
residence time of the fluid particles in the chaotic region. 

The unsteady flows over complex geometries that we described in this paper show the same basic 
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mechanisms of mixing as for the OVP flow. Besides the fluid permanently entrained by the vortex, 
and the fluid flowing unperturbed downstream, the scenario is enriched by fluid particles temporarily 
trapped by the vortex and then detrained, before flowing downstream. This phenomenon, known as 
transient chaos, was acknowledged and investigated a few years ago by Pentek, Tel & Toroczkai (1995) 
in the context of the advection problem of passive tracers in the velocity field of vortex pairs. 

We focus on the phenomenon of heteroclinic chaos as a large scale mechanism of intermittent release, 
causing intense concentration fluctuation of passive tracers. The entangling of the heteroclinic orbits is 
a consequence of the unsteadiness of the recirculating region and, in particular, of the oscillatory motion 
of the reattachment location. This oscillation has been observed in large eddy simulations as well as in 
two and three-dimensional direct numerical simulation of high Reynolds number flows over a backward- 
facing step (Driver, Seegmiller & Marvin 1983, 1987; Le, Moin & Kim 1997 and references therein). The 
periodic movement of the reattachment location is caused by the formation and detachment from the 
step of large-scale coherent vortical structures. The discrete vortex method that we applied to the snow 
cornice reproduces the same phenomena: figure 3a shows the fluctuations of the total circulation during 
the transient, due to the periodic detachment of clusters of vortices from the recirculating region. The 
phenomenon is evident during the transient because no perturbations were introduced in this simulation, 
but a natural oscillation of the centre of vorticity takes place during the transient, as shown in figure 3b. 
The transient is characterized by the roll-up of the vortex sheet which is shed into the fluid (see figure 2), 
consistently with other discrete vortex models of reattached flows (Cortelezzi, Leonard & Doyle 1994) 
and DNS of the backward-facing step flow (Le, Moin & Kim 1997). 

As shown below, the unsteady single vortex model is able to display the periodic detachment of 
coherent structures, even though these structures are not vortices but fluid particles trapped into the 
lobes of perturbed heteroclinic orbits. 

We present the results of two simulations showing the intermittent release of tracers. The first 



refers to the snow cornice generated by the modified Ringleb transformation (19). We perturbed the 
flow by displacing the vortex by e = 0.005 and letting it move freely along its periodic orbit. Passive 
tracers are released at a constant rate from a point upstream of the edge and close to the wall; their 
time evolution is described by integrating the governing equations ( |4l| ) by means of a fourth-order 
Runge-Kutta algorithm. 
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The time evolution of the plume is obtained by recording the position of the tracers at fixed time 
intervals. Also, we record the scalar concentration, averaged over a small control volume, as a function of 
time. Figure 8 shows the distribution of scalars in the field and their concentration in a sampling volume, 
marked by a shaded rectangular area, at two different times. The sampling volume is downstream. The 
periodic character of the concentration is demonstrated by its time history, plotted in the upper portion 
of each frame. The wave form of the concentration was found to be dependent on the location of the 
monitoring point. 

A second example of intermittent release is given by the simulation performed for the geometry ob- 
tained by the modified Schwartz-Christoffel mapping (^), i.e. the backward- facing step with smoothed 
corner. The unsteadiness of the flow is simulated by imposing a random perturbation with zero mean on 
the position of the vortex in equilibrium. Again, the concentration at a point downstream as a function 
of time in figure 9 illustrates the intermittent formation of lobes with high density of tracers which are 
responsible for the concentration fluctuations. Because of the random nature of the perturbation, the 
concentration is not periodic in time. 

The strongly evident feature of both the simulations is that the perturbation affects the detachment 
region slightly but the reattachment region quite strongly: the upstream stagnation point moves imper- 
ceptibly while the downstream stagnation point has a patently larger movement. This is consistent with 
the behaviour of real flows around edges, where viscosity effects constrain the separation to be located 
at the edge, making the separation region rather insensitive to perturbations. From the analytical point 
of view the two stagnation points in the (^-plane have oscillations of comparable amplitude, however the 
mapping onto the physical z-plane alters the displacements. The lengths in the z-plane are proportional, 
to first order, to the absolute value of the mapping derivative \dz/dC,\, whose minimum value is at the 
edge (for a sharp convex edge dz/dC, — 0). This analytical property provides a reasonable extension of 
the above mechanism of scalar intermittent release to the geometries characterized by high curvature 
edges, as well as by separation and reattachment of the flow. 
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5 Conclusions 



This study has examined the effects of the unsteadiness of reattached flows on the large-scale dynamics of 
dispersion of passive scalars. We considered the two cases of flow past a snow cornice and a backward- 
facing step, and we used point-vortex models to simulate unsteady separated shear flows. A new 
technique has been described to determine the streamfunction in a fleld of constant vorticity. 

A rotational flow over the above geometries reproduces the separation, even though the location of 
the separation point must be corrected by adding a point vortex to the flow, in order to have separation 
at the point of maximum curvature. The classical Ringleb mapping for snow cornices, and a Schwartz- 
Christofli'el mapping for step flows have been modified by smoothing the corners. Such modification 
is necessary to introduce the unsteadiness of the recirculating bubble, in that the separation point is 
allowed to move around the point of maximum curvature, while having the flow still satisfy an unsteady 
Kutta condition. 

The simulation performed for the snow cornice by a more sophisticated multiple-vortex model de- 
tected the oscillation of the center of vorticity, thus showing that the oscillatory motion of the reattach- 
ment point is associated with the unsteadiness of the recirculating bubble. The same kind of unsteadiness 
was imposed in the simulation by the simpler single-vortex model, by perturbing the equilibrium po- 
sition of the vortex. As the multiple-vortex simulation shows, the flow is characterized by a periodic 
release of clusters of vortices from the recirculating region behind the crest, in agreement with the direct 
numerical simulations of the backward- facing step of Le, Moin & Kim (1997). The same mechanism of 
releasing vorticity holds for passive tracers. 

Two different perturbations were applied to the system: a periodic, natural perturbation obtained 
by displacing the vortex from its equilibrium location, and a random motion with zero mean imposed on 
the vortex in equilibrium. The scalar concentration averaged over small control volumes was calculated 
as a function of time. For the case of periodic vortex motion, our model detects the intermittent release 
of blobs of tracers with the same period as the movement of the reattachment location (i.e. the same 
period as the vortex motion). For the case of random oscillation of the recirculating bubble, the release 
of coherent structures of tracers was still observed, although the wave form of the concentration is not 
longer periodic. 

Intermittent dispersion is associated with the phase-space lobe dynamics (Rom-Kedar, Leonard & 
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Wiggins 1990, Wiggins 1992), which is an intrinsically intermittent phenomenon. Since very small 
perturbations arc sulficicnt to start the above large-scale phenomenon, we expect intense concentration 
fluctuations to appear downstream whenever a flow undergoes separation and reattachment. 

Much of this work was done while P. Pranzese was a postdoctoral fellow at CSIRO Atmospheric 
Researh at Aspendale, Australia. We express our appreciation to Dr. Mark Hibberd for his thorough 
comments, as well as for his continuous support. We are also grateful to Dr. Mike Borgas for many 
helpful discussions and for his encouragements. 



Appendix A 

The Hamiltonian of a vortex 

The Green function G(C) ^o, ??o) for a flow region TZ. confined by a line C that extends to infinity is 
defined as 

G(C,r?;Co,%) = 9{i, V, Co, Vo) + ^ log v'($ - ?o)^ + {v - mf 

where ^(C, Co, %) is such that ry;^o,%) = on C, and is harmonic with respect to (C, ??) on TZ. 
Moreover, the reciprocity property of Green functions, i.e. G(C, Co, %) = (j(^o,??o;C,^), leads to the 
relationships: 

d d 

■^g{^o,m;^o,m) = 2 liin ^5(C,??;Co,%) 

d d 
;^5'(Co,%;Co,%) = 2 lim —5(^,77; ^0,%) • 
dr]o C^Co or] 

The streamfunction ■^i'(C, Co, ??o) due to a free vortex located in ^ = (^o must satisfy the boundary 
condition ij) = const on C, and has to be harmonic in the whole region except at Co, where it behaves as 
(7/277) log -\/(C — Co)^ + (*? ~ ^o)^- Therefore, the streamfunction is related to the Green function by: 

V'(C,^;Co,%) =V'i(C,^)+7G(C,r?;Co,%) (46) 

where V'l is the streamfunction of other flows possibly superimposed on the vortex. 

A vortex does not induce velocity on itself. Therefore, it moves as a particle in the flow fleld from 
which the vortex induction has been removed. We can calculate the complex conjugate velocity of a 
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vortex as: 



where w is the complex potentiaL According to equations (p5|), the velocity of a vortex can be expressed 
as: 

• dH' dH' 
dm) d£,o 

where the Hamiltonian H' is given by: 

(49) 



H' = V^i + . 

For instance, in our case the curve C and the flow region TZ are the real axis and the positive imaginary 
half-plane respectively, and the flow is the superimposition of a uniform flow and a vortex. Therefore 
the vortex Hamiltonian H' (^) becomes: 



An 



(50) 



The Routh rule 



When the motion of a vortex is conformally mapped onto the physical z-plane, its Hamiltonian has to 
be corrected according to the so-called Routh rule (Routh 1881). The following is a brief account of its 
derivation. 

The complex conjugate velocity of a vortex in the z-plane is 



■* ■ ■ ■ r I'dw 7 1 



zo \ dz 2-Ki z — zq 



This can be written as: 



in = lim 



dw 7 1 \ d( 7 



that is 



d(^ 2ni C — Co / dz 2ni 



1 



iC- Co) 



dz 

dC 



z - zo 



^ ^^dCo _ 7 dCo ^ ^^dCo _ _7 d / dzg 

° ° dzo 2m ^ f dzo\^ ° dzo Am dzo \ dCo 
dC'o 



(51) 



(52) 



(53) 
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Recalling equations (Eq), it follows that 



dH' d£,a dH' drjo .fdH'd£,o dH' dija 



d£,o dyo dfjQ dyo 



d 

dyo 



dxo^ 



V d^o dxo 
dzo 



drio dxo 



— f— log 

Ait \dyo 



dzo 



dCo 



z-— log 

dxo 



dzn 



T 1 



and finally 



H = H' 



J_ 

in 



log 



dzQ 



(54) 



(55) 



The relationship ( [49[ ) between a vortex Hamiltonian and a Green function for a bounded simply con- 
nected region was first given by Masotti (1931) (but see also Caldonazzo 1931 and Pclosi 1926). Succes- 
sively, Lin (1941) extended the Masotti theory to multiply connected domains containing more than one 
vortex. Surprisingly, Ringleb (1961) seems to have been unaware of the previous studies, as he deduced 
independently the Routh rule. However, he considered that it was only applicable for mappings based 
on many- valued functions. 
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FIGURE CAPTIONS 



Fig. 1: (a) The snow cornice obtained by the Ringleb mapping; (b) the backward- facing step 
obtained by the Schwartz-Christoffel mapping. 

Fig. 2: Simulation of the flow past a snow cornice by a multiple vortex method. The solid circles 
represent the vortices shed throughout the transient. The vortices strengths are proportional to the 
symbol sizes. 

Fig. 3: (a) Total circulation for the multiple vortex method (solid line) and circulation for the single 
vortex model (dotted line), as a function of time; (b) distance of the centroid for the multiple vortex 
method from the equilibrium point of a single vortex. 

Fig. 4: Streamlines and velocity vectors for a separated shear flow for the blunt snow cornice. 

Fig. 5: Pattern of the possible trajectories of the vortex. The equilibrium point is marked with a 
square. 

Fig. 6: (a) The separatrix streamline connecting two fixed points (P~,P+) in the unperturbed 
case; (b) as a consequence of the perturbation, the streamline splits in two branches which intersect 
each other: unstable manifold (bold line) and stable manifold (thin line). 

Fig. 7: Streamline pattern of a vortex pair flow in the frame of reference moving with the vortices. 

Fig. 8: Distribution of scalars in the field behind a snow cornice and concentration time series 
sampled in the volume marked by the shaded rectangular area downstream, at two different times. The 
concentration time series is plotted in the upper portion of each frame. The system is periodically 
perturbed by the vortex autonomously oscillating around its equilibrium location. 

Fig. 9: Distribution of scalars in the field behind a backward facing step with smoothed corner and 
concentration time series sampled in the volume marked by the shaded rectangular area downstream, 
at four different times. The concentration time series is plotted in the upper portion of each frame. A 
random perturbation is imposed on the vortex equilibrium location. 
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